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PREFACE 


This  report  is  intended  to  be  the  first  of  a  series  that  will  at¬ 
tempt  to  incorporate  some  of  the  latest  theoretical  techniques  into  a 
unified  but  concise  treatment  of  crystal  field  theory  and  the  analyti¬ 
cal  methods  used  therein.  The  series  of  papers  assumes  familiarity 
with  the  fundamentals  of  quantum  mechanics,  and  a  firm  background  in 
general  physics  and  mathematics. 

The  book  of  E.  U.  Condon  and  G.  H.  Shortley,  "The  Theory  of  Atomic 
Spectra,"  first  published  in  1935  and  reprinted  as  late  as  1959,  pro¬ 
vides  a  solid  basic  reference  to  the  field,  Since  the  time  it  was 
first  published,  there  has  been  an  increased  use  of  group  theoretic 
methods  (Bethe  used  these  methods  in  a  paper  as  early  as  1929)  and 
"Racah  Algebra"  in  the  analysis  of  the  atomic  spectra  and  their  modi¬ 
fication  due  to  crystalline  fields.  The  book  by  M.  E.  Rose,  "Elemen¬ 
tary  Theory  of  Angular  Momentum,"  provides  an  excellent  account  of  the 
Clebsch-Gordan  (vector  coupling)  coefficients  and  the  Racah  (recoupling) 
coefficients,  the  manipulation  of  which  constitutes  the  so-called  "Racah 
Algebra." 

The  standard  reference  to  group  theory  is  E.  P.  Wigner's  book  by 
that  name,  but  there  exist  references  that  are  more  readable.  Examples 
of  these  are  J.  S.  Lomont,  "Application  of  Finite  Groups,"  and  M.  Hammer- 
mesh,  "Group  Theory." 

Although  these  latest  references  are  readily  available,  there  are 
several  examples  in  the  literature  where  the  methods  employed  are  not 
quite  up  to  date.  There  is  failure,  in  some  instances,  to  recognize  the 
unified  treatment  of  quantities  which  have  been  given  various  symbols  by 
different  authors  but  are  simply  related  to  C-G  and  Racah  coefficients. 

As  a  result,  some  unnecessarily  involved  and  lengthy  expressions  have  re¬ 
sulted.  Now  with  the  availability  of  complete  tables  of  3-j  and  6-j  sym¬ 
bols  which  are  proportional  to  the  Clebsch-Gordon  and  Racah  coefficients 
respectively,  the  expression  of  results  entirely  in  these  quantities  is 
at  least  desirable,  if  not  necessary. 
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ABSTRACT 


This  paper  is  intended  to  be  the  first  of  a  series  that  will  at¬ 
tempt  to  incorporate  some  of  the  latest  theoretical  techniques  into  a 
unified  but  concise  treatment  of  crystal  field  theory  and  the  analyt¬ 
ical  methods  used  therein.  The  emphasis  in  this  paper  is  on  angular 
momentum  states  and  the  so-called  Clebsch-Gordan  and  Racah  coefficients 
which  arise  in  the  handling  of  such  states. 

In  order  to  illustrate  the  use  of  the  formalism  introduced  in  con¬ 
nection  with  the  angular  momentum  states,  a  simple  problem  in  quantum 
mechanics  is  solved  in  detail.  The  problem  brings  out  some  features 
of  atomic  systems  such  as  the  Zeeman  effect  and  spin-orbit  coupling,  and, 
in  addition,  it  illustrates  the  use  of  first  and  second-order  degenerate 
perturbation  theory. 

1 .  INTRODUCTION 


In  the  treatment  of  quantum  mechanical  systems  with  bound  states  due 
to  central  potentials,  it  is  convenient  to  describe  the  system  in  terms 
of  angular  momentum  states  (ref  1).  Even  in  cases  where  the  spherical 
symmetry  is  destroyed  by  the  imposition,  say,  of  an  external  magnetic  or 
electric  field,  angular  momentum  states  in  many  cases,  still  serve  as  a 
useful  basis  in  which  to  represent  the  system.  For  this  reason,  the 
earlier  sections  of  this  paper  describe  the  essential  features  of  angular- 
momentum  states  and  their  associated  properties.  The  later  sections  of 
this  paper  utilize  some  of  these  properties  in  the  analysis  of  an  illustra¬ 
tive  problem  and  introduce,  in  addition,  the  techniques  of  elementary  per¬ 
turbation  theory. 

2.  ROTATIONS  OF  ANGULAR  MOMENTUM  STATES 


Let  us  denote  an  angular  momentum  state  by  |jm\.  The  arguments  as¬ 
sociated  with  this  state  indicate  that  it  is  an  eigenstate  of  JB  and  J 
thus,*  Z 

J*  |jm^  =  J(j+1)  | Jm ^  (1) 

Jz  =  m  (2) 


where  «fa  and  Jz  are,  respectively,  the  square  of  the  total  angular  momen¬ 
tum  operator  and  its  z  component. 

The  total  angular  momentum  operator  may  be  defined  in  terms  of  the 
three  dimensional  rotation  operator  R  (^,9), (ref  2 )'  where  ^  is  the  unit 


*  The  units  *tf  =  c  =  1  are  used  throughout  this  report. 
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Vector  about  which  the  rotation  through  an  angle  0  is  effected.  An 
alternative  description  of  the  rotation  is  provided  by  the  Euler  angles 
(ref  3)  so  that  R  (ft,0)  may,  equivalently,  be  denoted  R  (ap  if).  In 
terms  of  such  a  rotation  operator,  the  angular  momentum  operator  is 
defined  by 


R(n,9)  |  jm  ^  =  e  ^  ^  |jm^  (3) 

or,  equivalently, 

R(apY)  |jm^  =  e‘iJza  e"iJyP  e~iJ*Y  | jm  (4) 

Assuming  the  |jm^  to  have  a  normalization,* 

<j'm'|W>  =  6JJ#  (5) 

the  matrix  representation  of  R  (a  p  V)  in  that  basis  is  defined  thus, 

/j'm'  |R(a  p  V)  |jm^>  =  Sjj/  D^/m  (a  0  Y>  (6) 

The  rotation  operator  is  diagonal,  in  j  since  the  angular  momentum  of  a 
system  is  not  altered  if  the  system  is  viewed  from  a  rotated  frame. 

The  projections  of  j,  however,  are  obviously  not  conserved  under  such 
an  arbitrary  rotation.  Using  the  completeness  relation  for  the  states 
|jm>  ,  i.e. 

lJln^  ^Jml  =  i  (7) 

Jm 


one  may  determine  the  result  of  rotating  an  arbitrary  state.  The  re¬ 
sult,  using  (6)  and  (7)  is, 


R(a  P  Y) 


1  •  Wap  Y) 


£  jj'm^  /m/  | R(OC  0  Y)  |jm^> 


J  * 


- 1  »i'.  p  v> 

_  /  r 


(8) 


*  For  those  not  familiar  with  the  .Dirac  bra-ket  notation,  associate  wave 
functions  with  the  |jin^  and  dOfj’^/  Pf  jm  with  'm' |p|  jm^  for 

the  matrix  element  of  an  arbitrary  operator  P  .  Further  information 
may  be  obtained  in  "Quantum  Mechanics,"  E.  Merzbacher,  p.  306,  (ref  1) 
or  "The  Principles  of  ^lantum  Mechanics,"  P.  A.  M.  Dirac  (ref  10). 
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Thus  one  has  expanded,  into  the  set  of  states  that  form  the  basis,  the 
result  of  rotating  the  state  | Jm  ^  (which  itself  is  one  of  the  '’members" 
of  the  basis)  through  angles  defined  by  a,  (3,  Y.  To  avoid  confusion,  it 
is  here  emphasized  that  the  operator  R  (a  0  Y)  as  defined  above  operates 
on  the  physical  system  described  by  the  ket  and  physically  rotates  that 
**  system.  After  such  a  rotation,  the  physical  system  is  not  (in  general) 

in  a  pure  state  insofar  as  the  basis  of  states  (the  | j m  ^> )  is  concerned. 
It  is  in  a  state  of  impure  z  projection  as  evidenced  by  its  description 
w  as  a  superposition  of  pure  states  with  various  m  values  as  given  by  (8). 

In  many  places  in  the  literature  one  will  read  that  the  coordinate  system 
is  being  rotated  and  in  many  instances,  it  is  not  clear,  to  which  coor¬ 
dinate  system  reference  is  being  made.  The  interpretation  given  above 
seems  to  be  less  confusing. 

It  is  convenient  at  times  t6  write  the  explicit  dependence  of 

/  (a  (3  Y)  on  its  arguments  a  and  Y.  That  this  is  possible  is  easily 
m  in 

seen  if  the  expression  (4)  is  substituted  into  (6).  The  result  is 


D^'b^oPY)  =  j  e  iJ*a  e  iJy^  e  |jm^ 

(9) 

.  <V|  .-‘¥ 

since  the  Jz  can  be  replaced  by  the  eigenvalues  of  the  states  on  which 
they  operate.* 

The  matrix  element  in  (9)  is  denoted  by  the  symbol  d^ /  (0),  so  that 
.  mm 


Di'n(apy)  ■  e”im/a  «£/m(0)  e~imy 

in  la  in  m 

(10) 

where 

di/m(P)  E  <(jm'ie’1JyP 

(10a) 

To  fix  some  of  these  results  and  to  relate  the  to  known  functions, 
consider  the  rotation  of  a  system  described  by  a  Legendre  Polynomial 
P,(cos  9Q),  (ref  1,  4).  Now,  0  can  be  considered  as  the  angle  between 
two  unit  vectors.  It  is  simply  a  function  of  the  dot  product  of  these 
two  unit  vectors  which  are  purely  intrinsic  to  the  system,  and,  thus. 
Independent  of  any  external  reference  frame.  Therefore,  one  must  have, 


*For  example,  such  a  requirement  can  be  thought  to  define  what  is  meant 
by  e  iJza,  viz.,  e  iJ*a  |jm)  ■  e  ima|jm^if  |jm^>  =  m  |jm^  • 
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R(C$Y)  p.  (COS  0  )  =  P.(cos  0  )  (11) 

t  o  C  o 

Independent  of  the  rotation,  From  (A15)  Appendix  Aj,  P^  may  be  expressed 
in  terms  of  the  spherical  harmonics,  Y^  (ref  4),  thus  enabling  a  differ¬ 
ent  interpretation  of  RP  .  Since  the  Y^m  are  angular  momentum  states, 
we  may  apply  formula  (8)^0  obtain 

R(ctfJV)  Pt(co,  »o)  .  R(OPV)  Yw(So,0) 


(COY)  Yim,(0,cp) 


(12) 


In  the  above,  the  Y.„(8  ,0)  refers  to  a  particular  orientation  of  the  phys- 
7  <0  o 

ical  system  with  one  of  the  unit  vectors  in  the  x-z  plane  (hence  -  0) , 
and  the  other  specifically  along  the  z  axis,  so  that  the  angle  between  the 
two  is  the  0O  as  required.  The  angles  0,  cp  are  the  spherical  coordinates 
of  the  former  unit  vector  after  rotation.  Using  the  invariance  of  P^ 

(cos  0  )  as  expressed  by  (11)  and ' comparing  the  result  (12)  with  the  addi- 
tional°formula  for  the  spherical  harmonics  given  by  (A15),  one  concludes 
that 

to™’  -(sTri)1*  Y«  <e’0)  • 

That  (O0Y)  is  independent  of  Y  4s  seen  by  (10).  The  expression 
mO 

is  an  important  relationship  and  is  used  later. 

The  D^  are  thus  shown  to  be  simply  related,  for  certain  values  of 
their  arguments,  to  the  more  familiar  spherical  harmonics.  Some  of  the 
properties  of  the  I)J  are  listed  in  Appendix  C  for  convenience. 

3.  ADDITION  OF  TWO  ANGULAR  MOMENTUM  STATES 


(13) 

(13) 


When  adding  two  angular  momentum  states  to  obtain  a  combined  state 
that  has  a  definite  total  angular  momentum  (ref  2),  one  necessarily  loses 
some  information  in  the  process.  Initially,  one  possesses  a  "product" 
state 

IVl^  lJ2m2^>"  lJimiJ2m2^>  (14 

which  is  defined  by  four  quantum  numbers.  It  is  desired  to  obtain  (by 
some  superposition  of  such  states)  a  resultant  state  ^  that 

has  a  well-defined  total  angular  momentum  and  projection  as  indicated, 
and  in  which  the  individual  angular  momenta  are  preserved  as  good  quan¬ 
tum  numbers.  The  transformation  can  formally  be  effected  thus. 
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(15) 


^"^2-2 


The  matrix  element  is  proportional  to  6.  ./  so 

J2J2 


that  the  final  form 


lJMJlJ2^>  =  ^  ^  jlmlJ2m2  I  JMJ1J2^  lJimiJ2m2/> 


The  coefficients  given  by  the  matrix  element  above  are  generally  called 
Clebsch-Gordan  (C-G),  vector  coupling,  or  Wigner  coefficients  and  con¬ 
stitute  an  orthogonal  transformation.  Frorti  here  on,  the  notation  to  be 
used  for  them  is 


^j1<®1)  W  lJ(M>^  8  <^JimiJ2m2  lJMJlJ2^>  (17) 

denoting  that  angular  mementa  j^  and  j^  with  projections  and  are 

added  to  produce  a  resultant  total  angular  momentum  J  with  projection  M. 
Other  useful  properties  of  these  coefficients  are  listed  in  Appendix  B. 

In  the  simplest  one,  the  coefficient  vanishes  unless  m^  +  =  M. 

Therefore,  the  sum  in  (16)  over  one  of  the  m^  can  be  performed  to  obtain, 
in  the  new  notation, 

lJMJ1J2/>  =  I<\J1<m)  j2(M',n)  lJ(M>  y  (18) 

m 

In  this  form,  the  information  which  has  been  "lost"  in  order  to  gain 
information  about  the  total  angular  momentum  is  made  apparent.  In  summing 
over  m,  we  have  lost  information  regarding  the  individual  projections  m 
and  mg  written  m  and  M-m,  respectively.  In  other  words,  in  order  to  ob¬ 
tain  a  state  with  well  defined  total  angular  momentum  J  one  must  superim¬ 
pose  states  of  all  possible  projections  m^  and  m consistent  with  the  re¬ 
quirement  m  +m„  =  M.  The  amplitude  with  which  each  such  combination  enters 

into  the  superposition  is  given  by  the  Clebsch-Gordan  coefficient  as  shown 
in  (18) . 

An  immediate  application  of  (18)  can  be  made  to  the  electromagnetic 
field  of  any  such  vector  field  (ref  4).  The  orbital  angular  momentum  in 
such  a  field  is  given  by  the  Y,  apd  the  intrinsic  spin  of  that  field  is 
one,  hence  it  is  described  by  the  5?  ,  the  spin-one  spinors  (ref  4,  5). 
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i 
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The  Y.  Is  an  angular  momentum  state  of  angular  momentum  and  projection 
m.  Tne  X  is  an  angular  momentum  state  of  angular  momentum  one  and  pro¬ 
jection  Hence,  according  to  (18) 


U>,)  \ 

(i 


(19) 


where  the  equation  defines  the  vector  spherical  harmonics, 


These 


are^  wave  functions  of  total  angular  momentum  J  and  projection  M,  and  are 
useful  to  describe  the  radiation  from  atomic  states  that  are  expressed  in 
terms  of  total  angular  momentum  states. 


Another  example  that  is,  perhaps,  more  familiar  is  the  coupling  of 
two  Pauli  spinor  (spin  pne-half)  ^tates.  Let  represent  a  Pauli 
spinor  (ref  6)  where  p  =  ±  1/2  £hen, 


J 


u 


I  «•*>>  >W  V 


and 


♦oo  =  5  l°<0>>V  V 


(20) 


(21) 


The  states  given  by  (20)  are  the  familiar  triplet  states  for  two  parti¬ 
cles  usually  denoted  X*,  and  the  state  given  by  (21)  is  the  singlet  two- 

particle  state,  usually  denoted  xS • 

The  addition  formula  for  two  is  obtained  by  a  consideration  of  (18). 
Rotate  the  physical  systems  related  by  (18)  with  the  operator  R(oPY)-  "One 
may  write 


R(QJ3Y)  IJMj^jj^  =  J1(m)j2(M-m)  |j(M))r  R(QpV)  jj^^ 


'  RWV)  | J2M-in^>  ,  (22) 

since  the  state  |j.mj  M-m^  is  a  product  of  two  single  system  states  ac¬ 
cording  to  (14).  Taking  the  matrix  element  with  the  adjoint  equivalent 
of  (18)  and  referring  to  the  definition  of  I)j  given  by  (6),  one  obtains, 

(C^Y)  =  X<CJi<m)J2<ia_m) 

mm/ 

•  D^,  (apY)  (C0Y)  (23) 

mm  M  -m  M—m 
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The  expression  above  is  easily  inverted  i£  one  operates  on  both  sides 
with 


£  |J(M)^>  /j1(^/)J2(M/V)  |jy(M')^>  (24) 

Jj' 

and  uses  the  property  of  the  C-G  coefficients  given  by  (BB)  Appendix  Bi 
The  result  is 

|J(M)N  ^J1(|a/)j2(M/-m/)|  •KM/)^>dJm/  (af3Y) 

J  j,  J2 

•  V  <QPY>  DM-u  i'V<QPY>  (25> 

If,  now,  p'  and  m'  are  set  equal  to  zero,  the  expression  above  reduces 
(applying  (13)  and  taking  complex  conjugates)  to  the  expression  (A22) 
Appendix  A  involving  the  spherical  harmonics. 


The  result  (25)  enables  one  to  evaluate  Integrals  involving  three 
(or  more)  J>J  .  The  orthogonality  of  the  is  expressed  by 


(f>  do)  B3*  (O0Y)  D^/  /(QfSY) 
J  mu  m  |a 


Sir3  t 

2j  +  1  jj' 


(26) 


where  du)  s  da  d(cos  p)  dY.  Therefore,  using  (25),  it  is  simple  to 
deduce 


$ 


-C«*V)  °dx..  (a^Y)  °i2n  (QpY)  ■  57iT 

MM  m^  m2^2 


<CJl(mi)J2(m2) 


•<^J1(^1)j2(u2)  |j(M')^>  (27) 


where  the  Kronecker  deltas  involving  the  projections  are  implicit  in 
the  C-G  coefficients.  For  M/  *  =  \i„  =  0,  one  may  obtain  by  the  use 

of  (13),  the  fdllowing  expression  for  tne  integral  of  three  spherical 
harmonics, 


$ 


<m  Y 


1 M 


(9,W  Y.  „  (®>®>' 

Vl 


y.  _  (®,«P) 

l2m2 


(24  +i)(a*a+i]"B 

.  4V(2bf'l)  \41(®1)-f  2(m2)  |L(M)^ 


•<41(0)42(0)  |l(Q)^ 


(28) 


where  dfl  a  dtp  d(cos  9).  It  is  an  often  used  fact  in  the  theory  of 
atomic  spectra  that, 

^(0)^(0)  |L(0)^>  oe  £i  +  (-)>tl+'VL]  (29) 

i.e.,  it  vanishes  unless  41+-t2+L  is  an  even  integer. 
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In  the  manner  demonstrated  above ,  one  may  proceed  to  evaluate  the 
integrals  involving  an  arbitrary  number  of  or  spherical  harmonics. 

The  integrals  will,  of  course,  be  expressed  in  terms  of  the  C-G  coeffi¬ 
cients  . 

In  many  ways,  it  is  advisable  to  become  familiar  with  D3 .  They  are 
referred  to  as  the  irreducible  representations  of  the  three  dimensional 
rotation  group,  a  three-parameter  continuous  group.  Although  their  nor¬ 
malization  (26)  is  not  as  simple  as  that  of  the  spherical  harmonics,  in 
certain  applications  their  use  is  to  be  preferred.  Such  matters  will  be 
considered  in  a  later  paper. 

4.  ADDITION  OF  THREE  ANGULAR  MOMENTUM  STATES 


The  method  of  coupling  three  angular  momentum  states  (ref  2)  to  ob¬ 
tain  a  total  angular  momentum  state  proceeds  logically  from  expression 
(18).  The  procedure  is  to  couple  two  of  the  states  (any  two)  to  obtain 
a  resultant  two-particle  total  angular  momentum  state,  and  then  to  cou¬ 
ple  this  resultant  state  with  the  third  state,  again  according  to  (18). 

To  illustrate,  suppose  the  three  states  are  >  | J2m2  ^  and 

^3m3  ^  '  CouPle  the  first  two  to  obtain  ' 

| ,)m j i j 2^  =  Jl(ml)J2(m2)  :  lJ2m2/>  (30) 

remembering  that  the  sums  are  restricted  since  /i  (■. ) |j(m)/ 

■  6m1+m2,m-  Then’  COUple  *°  ^3*3}  *°  obtaln> 

|JRJJ3>  =  lJ3m3^  (31) 

mm3 

If  (30)  is  substituted  into  (31),  then  result  is 

|jM(J1J2)(jj3)>  =  |jCm)^><^j(m)j3(m3) 

mlm2mm3 

'  lJ3m3> 

Each  of  the  individual  angular  momenta  as  well  as  the  intermediate  j  are 
good  quantum  numbers  (since  they  have  not  been  summed  out)  and  are  there¬ 
fore  included  in  the  resultant  state.  Because  of  the  nature  of  the  C-G 
coefficients,  the  sum  over  projections  collapses  to  a  sum  over  only  two 
indices . 


l««M>^|Ji»1>|Jami^ 

(32) 
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(33) 


*S 


^  J1(M)J2(m“M)  ^J(m)J3(M-m)  |  j(M)^> 

mp 


The  order  in  which  the  individual  angular  momenta  appear  in 
IJMCJjjjKjJg)  y  is  very  Important  because  it  indicates  the  order  in 

which  the  states  were  coupled.  For  example.  I 'j. )\  de¬ 
fines  ,  231/. 

mV' 

I 

04) 

The  order  in  which  the  single  particle  states  appear  in  the  summation  is 
immaterial  since  they  do  not  interact  amongst  themselves.  What  is  im¬ 
portant,  however,  is  the  order  of  the  angular  momenta  in  the  C-G  coeffi¬ 
cients.  It  is  emphasized,  too,  that  the  projections  of  the  individual 
angular  momenta  are  the  same  in  the  C-G  coefficients  and  the  states  that 
they  couple. 


It  is  frequently  of  interest  to  express  a  three-particle  state  that 
haB  been  formed  by  one  method  of  coupling  in  terms  of  states  coupled  in 
a  different  manner.  This  expression  is  determined  by  obtaining  the  fol¬ 
lowing  matrix  element, 


'j')  IJMCJjjgXjjg)^  a  ^J1(M)J2(m-M)|j(m)^><^J(m)j3(M-m)|j(»l)^> 


m 

m'|i/ 


•  ^j2<MVj3<»'V)|jVm')^^J,(m')j'<M'-m')|j'<M,)^>fiJ  jf6j  j  ,6  ^ 


'll-mjm'-p' 


(35) 


using  the  normalizations  given  by  (5).  The  product  of  the  three  Kronecker 
deltas  in  (35)  involving  the  projections  are  equivalent  to  m' 

ao  that  the  sum  over  m'  and  p'  can  be  performed.  The  result  is 
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')(J -J^>  I  J«(J1J2)  <JJ3)>  - «,  J.  ®j  J '  «J,  -  >m. 

X  J-  Z  m  3  3 

^j2(n>-4)J3CM-m)  |j'(ll-p)  y  |j;(M)^ 


(36) 


In  order  to  allow  the  use  of  (B8)  Appendix  the  C-G  coefficients 
Involving  j'  and  J '  must  be  recoupled  so  that  they  assume  the  form  of 
the  first  two  coefficients.  To  effect  this,  the  so-called  Racah  coef¬ 
ficients  (ref  2),  W  (a  b  c  d;  el),  must  be  used.  These  coefficients 
accomplish  the  task  desired,  as  follows: 

^a(p)b(m-p)  |e(«)^  ^e(m)d(M-m)  |c(M)^  - 

=  ^  [(2e+l)(2f+l)]1/8  W(abcd;ef) 

•  ^b(m-p)dOHa)  |f(M-^)^>  ^a(p)f(M-p)  |c(M)^  (37) 

In  order  to  apply  (37)  to  our  problem,  first  j  and  J3  must  be  inter¬ 
changed.  Using,  whenever  necessary,  formula  ( B5 ),  we  obtain, 

<^J2(m-M)J3(M-m)  |j'(M-p)^  <^j/(M-4)j1(p)  |j'(ll)^> 

=  (_)J2+J3  J  <^j3(M-m)J?(m-p)  Ij'Ol-p)^  </j  '(M-fiX^p)  |y<M)^> 

=(_)J2+J3'J  £[2j'+l)(2f+l)]^8  W(J3J2J/j1;j'f) 

f 

'  J2(m-t-i)J1(n)  |*(m))>  <(j3(ll-m)f(m)  |j'(M)^> 

=(-)J  '■3l”J^[2j/+l)(2f+i)]3/8  W(j3j2J/J1;j,f) 
f 

*  <(j1(^J2(m-H>  |f(m)^  <^f(m)j3(M-m)  |j'(M)^  (38) 


If  the  result  (38)  is  utilized  in  (36),  then  the  sum  over  p  may  be  per¬ 
formed  to  yield  according  to  (B8).  This,  then,  enables  the  sum 

over  f  in  (38)  to  be  performed,  fixing  f=j .  Then  the  sum  over  m  in  (36) 

can  be  performed  to  yield  6TT/.  The  final  result  is, 

JJ 
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(39) 


<jVo£'Mj'j{>  |J»<J1J2)(JJ3)>  »  »jij;  4j2j^  •*j3J '  V  V 

.  (-)J'Jl"j/[(2J+l)(2j'+l)]lA  W(J3J2JJ1;J/J) 
from  which  we  deduce 

|jM(j1d2)(jJ3>^  =  £(-)J"Ji“J  '[(aj+imj'+i)]1*  w(j3J2Jj1;j/j) 

y 


•  |jM(J2J3)(j/j1)^ 


<4p) 


The  Kronecker  deltas  in  (39)  could  have  been  predicted  at  the 
start.  In  each  of  the  states,  the  total  angular  mpmentum  and  its  pro¬ 
jection,  and  each  of  the  individual  angular  momenta  are  good  quantum 
numbers  and,  what  is  more  essential,  are  independent  of  the  mode  of  coup¬ 
ling.  Thus,  since  the  only  way  in  which  the  two  three-particles  states 
differ  is  in  this  latter  property,  one  can  expect  and  understand  the 
orthogonality  in  those  discrete  variables  embodied  in  the  Kronecker  deltas. 
The  intermediate  angular  momenta  j  and  j '  however  are  not  on  the  same 
footing  since  the  former  is  the  result  of  j  and  J2  coupling  in  a  pre¬ 
scribed  manner,  and  the  latter  is  the  result  of  a  coupling  between  j2  and 
j  .  Formula  (4)  expresses  the  fact  that  |jM( jjjg) ( J jg)  )  i»  *  mixture 
or  states  | JIKJ2J3) (J  ' Jl)  )  with  various  values  for  the  intermediate  angu¬ 
lar  momentum. 

The  symmetry  properties  of  the  Racah  coefficients  are  given  by  Rose 
(ref  2t  p  226).  (  71)ey  are  related  to  the  Wigner  6-J  smybols  in  the  fol¬ 
lowing  manner,  .given  by  Edmonds  (ref  5), 

a  b  e 

d  c  f 

Further  information  regarding  the  6-j  symbols  is  available  in  Edmonds,  and 
a  complete  table  of  their  values  for  angular  momenta  less  than  or  equal  to 
8  is  given  by  Rotenberg,  et.  al.  (ref  7). 

5.  APPLICATION  OF  THEORY 


W(abcd;ef)  =  (-) 


a+b+c+d 


The  foregoing  general  discussion  shall  now  be  utilized  to  calculate 
the-energy  levels  of  a  rather  simple  problem.  The  problefi  is  th#  7 
consideration  of  a  single  electron  in  a  coulombic  potential  (ref  6)  with 
spin  orbit  coupling  and  a  magnetic  field.  The  specific  Hamiltonian  is 

H  =  Hq  +  Hx  +  H2  (42) 
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where 


oa 

H  =  —  +  V(r)  with  p  =  -1  V 
o  m 


Hg  =  X  L  •  S 

Although  the  problem  Is  quite  simple ,  Its  solution  has  the  advantage 
of  illustrating  many  general  features  common  to  more  complicated  cases. 

The  problem  is  treated  first  by  perturbation  theory  (ref  1,  6) 
using  states  that  are  diagonal  in  Hq  +  with  H2  being  the  small  pertur¬ 
bation'  and  carrying  the  calculation  through  to  second  order  in  Hg . 

The  procedure  is  then  repeated  using  states  that  are  diagonal  in  H0  and 
Hg  with  Hj  considered  small.  Finally  the  problem  is  solved  exactly, 
and  the  result  is  expanded  to  second  order  in  Hg  and  then  to  second  or¬ 
der  in  Hi  to  show  the  agreement  of  the  results  obtained  by  perturbation 
theory  in  each  case  with  the  exact  result. 

When  the  potential  is  coulombic,  the  energy  levels  of  Hq  depend 
only  on  the  principal  quantum  number  n  (ref  6).  The  eigenfunctions  of 
HQ  are  then  of  the  form  where  l  is  the  orbital  angular  momen¬ 

tum  quantum  number  and  m  its  projection,  and  is  a  Pauli  spinor  re¬ 
presenting  the  spin  of  the  electron.  Only  the  case  n  =  2  is  discussed, 
hence  £  =  1,  0  and  ■  =  1,  0,  -1  and  ^  =  ±  1/2.*  The  energies  given  by 
and  Hg  will  then  be  relative  to  the  energy  level  associated  with 
n  =  2  for  the  operator  Hq  above. 


6.  CASE  I:  SpALL  Hg 

In  this  case  we  consider  the  magnetic  energy  much  larger  than  the 
spin  orbit  energy  and  treat  the  latter  as  a  perturbation.  The  magnetic 
field  will  be  taken  to  pbint  in  the  negative  z  direction,  so  that 

Hi  -  Hb  H(Lz  +  2V  (44) 

where  u  is  the  Bohr  magneton  (ref  8).  The  states  described  by  the 
spherical  harmonics,  Y^m,  are  diagonal  in  L^,  i.e.,  L^Y^  =  mY£m  and 

the  Pauli  spinors,  are  diagonal  in  S^,  i.e.,  Sz  X£i  =  For  con¬ 

venience  we  will  replace  m  by  M  -p  in  the  Y^m  above.  Then 


+  M  V  =  ^  H(M+m)  (45) 


♦The  number  of  states  giving  the  same  energy  (degeneracy)  for  an  electron 
in  a  coulombic  potential  is  2n8  (equal  to  8,  in  this  case). 


and  the  energy  due  to  the  magnetic  field  is 


E 

Mp 

=  (3(M+n) 

where  3  * 

then  obtain  the  i 

energy  levels  listed  in  the  Column  headed  by 

table  I . 

TABLE  I. 

ENERGY 

LEVELS  TO  SECOND 

ORDER 

IN  H2 

1st 

2nd 

1 

M 

M 

eM4 

order 

order 

1 

3/2 

1/2 

2^ 

X 

2 

0 

X 

X8 

1 

1/2 

-1/2 

0 

*2 

“23 

X8 

1 

1/2 

1/2 

3 

0 

23 

X8 

1 

-1/2 

-1/2 

-3 

0 

“23 

X 

X8 

1 

-1/2 

1/2 

0 

"2 

23 

1 

-3/2 

-1/2 

-2p 

X 

2 

0 

0 

1/2 

-1/2 

3 

0 

0 

0 

-1/2 

-1/2 

-3 

0 

0 

see  that 

ther?  are  eight 

states,  as  there 

should 

be,  with  3  1 

(46) 


MU 


the  energy  levels  two-fold  degenerate. 
The  perturbing  Hamiltonian  is 


XL-  S  =  XS  L  +  -  [L+  s”  +  L_  S+] 
z  z  « 

whe  re 

IT  =  L  ±  iL  ,  S*  *  S  ±  iS 
x  y  7  x  y 

We  will  need  the  relations  (ref  1) 

L±  vt«  -  W 


(47) 

(48) 

(49) 
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0 


(50) 


s“x  =  X  s+x  =  X  s'x  a  s+X 
p  -p  .  -p  H  -p  +p 

The  latter  formulae  are  special  cases  of  the  former  with  l  ■  1/2  and 
■  ■  li.  In  perturbation  theory  the  energy  of  a  given  state  m  to  second 
order  is  given  by  (ref  6) 

r  ,(S)  ("») 

M  1 1 v-  T-.  '  <“> 

N  ran 


wheje  ■  (j)  d  f)  $n  and  ^  are  eigenfunctions  of  the  un¬ 

perturbed  Hamiltonian.  Now 

x  “  '  t  TtM-n’V  ■  x  WM  *  5  <lV  *  L_8*>  Y«.-A  <M> 


•IM-p.  p 


-  X  <■-•*»>*  *  \  <LV*  (53) 

The  last  part  is  easier  to  carry  through  by  choosing  specific  values  for 
l,  M,  and  p.  For  =  1,  M  =.  3/2  and  p  i  1/2 

(bV  +  sV)  YnX1/3  -  o  (54) 

aince  Li^  ■  0,  and  ■  0.  The  same  result  is  obtained  for  t  ■  1, 

11  «  -3/2  and  p  a  -1/2  or 


X  8  '  L  Yl±l  X±l/3 


X  (2}  Yl±l  X±l/2 


Using  the  formulas  given  above ,  the  results  of  the  perturbation  acting 
on  the  remaining  state  Is 

X  8  •  L  YnX_1/a  ■  h3y11X_1/2  -  “  2  YHX-l/2  +  2  ^  Y10Xl/2 

H2Yl-lXl/2  "  "  2  Yl-lXl/2  +  2  ^  Y10X-l/2 

Vl0Xl/2  -^YllX-l/2  <M 

H2Y10X-l/2  *  2  ^  Yl-lXl/2 

V00Xl/2  -° 

H2Y00X-l/2  *  0 
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We  see  that  the  perturbation  in  no  case  connects  degenerate  states, 
so  that  the  second  order  energy  cap  be  calculated  by  the  simple  method. 
In  general  this  is  not  the  case  when  one  has  degenerate  states,  and  the 
problem  can  become  rather  involved  (ref  6,  p  155  ff.)  The  results  of 
the  calculations  for  first  and  second  order  are  shown  in  table  1. 


7.  CASE  II:  Hi  SMALL 

Ih  this  part  we  use  eigenfunctions  of 
call  these  functions  2^^.  Then  we  have 

-  ‘fA 

s‘fA  - 

The  reason  for  using  these  eigenfunctions  is  that  can  be  written 

Hg  a  X  L  •  S  «  5  [J*  -  La  -  Sa]  (58) 

j  s 

diagonal  in  the  functions  .  Further  since  Hq  is  degenerate 

angular  momentum  we  can  also  use  these  functions  as  eigenfunctions 
i.e. , 

»o  - «.  vyS  <59> 

This  choice  of  eigenfunctions  is  desirable  since  the  Hamiltoniap  Hq  +  Hg 

is  diagonal  in  that  basis.  In  case  I,  the  eigenfunctions  chosen  were 

Y.  X  ,  because  H  +  H,  (there  considered  as  the  unperturbed  Hamiltonian) 
-Cm  p7  o  1 

is  diagonal  in  that  basis.  Thus,  from  (57)  and  (58) 


and  is 


in  the 
of  H  , 


L,  S,  J  •  L  +  8  and  J  and 

z 


(57) 


X  L 


S J^JM  =  I  CJW+1)  '  ^+1>  -  WS+UtyS 


=  |  [j  (j+i)  -  ia+i >  -  (60) 

.slnde  S  =  1/2.  The  energy  levels  for  different  M,  and  J  values  are 
given  in  table  II,  in  the  column  headed  by  E  . 
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TABLE.  II. 

ENERGY 

LEVELS 

TO  SECOND  ORDER 

IN  % 

l 

1 

M 

E 

1st 

order 

2nd 

order 

0 

1/2 

1/2 

0 

0 

0 

0 

1/2 

-1/2 

0 

-p 

0 

1 

3/2 

3/2 

X 

2 

20 

0 

1 

3/2 

1/2 

X 

2 

2g 

3 

4 

27  X 

1 

3/2 

-1/2 

X 

2 

,.2§ 

3 

_i  g! 

27  X 

1 

3/2 

-3/2 

h\ 

2 

"20 

0 

1 

1/2 

1/2 

-X 

§ 

3 

-ll  g! 

27  X 

1 

1/2 

-1/2 

-X 

_g 

3 

._i  g! 

27  X 

Again  there  are  eight  states,  aa  there  should  be,  but  in  this  case 
two  energy  levela  are  two-fold  degenerate,  and  one  is  four-fold  degen¬ 
erate. 

We  now  consider  the  perturbation  H^.  This  can  be  expressed  as 

-  |T  •  H  -  ^,H  [Lz  +  28^]  -  0  [Jz  +  Sj  (61) 

Since  J  la  diagonal  in  the  above  representation  we  have  only  to  con¬ 
sider  tfte  perturbation  as  0SZ. 

Analagous  to  the  expression  (19),  we  can  write  the  following  for 
spin  one-half 

%  -  !<>“-“>  5'“>  !““>>  Wm  <“> 

9  n 

8  i<*>  !>“>>•*  VA  <63) 

*  n 

We  can  uae  expression  (B9,  Appendix  B)  to  invert  (62)  so  that 

WV  ■  j<“>  !•'»>  'fi'u  <M> 

J 
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Since  the  values  of  n  are  ±1/2  this  can  be  expressed  as 
M  *  j  (-l)1^  Then 

Sz^JM  *  HZ  <-1>1/3"M<(^<l,-^>  |  J(M)  >  4(M-m)  i(p)|jr/(M)>ljJ/|I  (65) 

0  M' 

Using  the  rules  for  manipulation  of  the  Clebsch-Gordan  coefficients  (Ap¬ 
pendix  B),  we  have 

<^0H0  i(p)  |J(M)>  -  (-1)1/2+^i||§  §H>  k(M-M)^  (66) 


Hence 


t ill  5»1*> i(M)  IJ'OO^., 


vJ+o  -l 


-  2  Vw+T  ^J/57TT  W(J  J  J'  i  (,f)  </j(M)  f(0)  |j'(M)^ 


f  J'p 


From  table  1  in  Rose  (ref  2),  we  obtain 


i(n)|f(0)> 


2'  2 


hence 


n/2<^  5 HO  |K0)>  -  1 

£<±HOj(M>  |1(0)>  <(|(-p)  i(n)  1 1(0)^-  6fl 


Finally  we  obtain 


~j  1-mrnTi)  »u  f 1  ~  ko>  k'(»>>|j-.u 

-j 


The  final  result  for  the  total  perturbation  is 
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(71) 


The  explicit  value  of  (j  |jj  lor  valued  of  j'  and  J  of  interest 

can  be  calculated  by  using  the  results  given  by  Rose  (ref  2,  Appendix  1) 
and  are  given  in  the  tables  below. 


TABLE  HI.  VALUES  OVd^Aj'  ;j) 

1=1  1=0 


V 

J  \ 

1 

2 

3 

2 

1 

2 

2M 

0 

3 

2 

0 

0 

. 

The  results  can  be  written  in  compact  form  as 

/  .  bm  (2j+d  n.i  {2  k  0/1 

<M  ”  2/+1  /Jlt  3  ^  4,11  6m±  1/2/  j/|1 


(72) 


which  holds  for  =  0,  J  =  1/2  and  for  £  =  1,  J  =  1/2,  j'  =  3/2  and  i  = 

Ja3/2,'j'  a  1/2.  Again  the  perturbation  connects  no  degenerate  states, 
so  that  the  energy  to  second  order  is  simply  proport  ional  to  the 
square  of  the  off-diagonal  elements,  i.e.  i  ^  j  ia 

4  Ba  a  3  l1  *** l/2 

give  four  elements,  — 
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These 

_  '  j'  ■  m*  i  /  y.  » 

4  B" 

^  twice  and  -  twice  and  are  shown  in  the 


last  column  in  table  II, 


EXACT  CALCULATIONS 


To  find  the  exact  diagonal  representation  we  can  reduce  the  calcu¬ 
lations  somewhat  by  observing  that  the  z  component  of  J  commutes  with 
the  total  Hamiltonian.  To  show  this  we  have 

H  ■  M  *  H  4  X  L  •  S*P  (L  +2S  )  +  XL  S  +  £  (L+s"  L~S+)  (73) 

E  Z  E  Z  a 

Where  L*  =  L  ±  1L 

*  y 

and  L  x  L  i  1L 

Then  [L^L1]  =  ±  L* 

and  [Sz,  S*]  »  ±  8± 
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Hence 


[(Lz+Sz),H  ]  =  |  [L+S~  «  L_S+  -  L+S"  +  jfs+]  =  0  (74) 

or  [J  ,  H]  =  0.  so  that  If  is  a  good  quantum  number.  To  utilize  this 
z 

fact,  we  use  the  functions  Y  X  in  place  of  Y,  X  .  Since  the  *’ 

t,ll-p  p  (n  fi 

Y  X  form  a  complete  set  of  functions,  and  since  J  commutes  with 

<44“  Z 

H,  we  can  form  the  eigenfunctions  of  H  from  a  superposition  of  the 

y*m-pV  that  18 

*  - 1  \  Vu  V  <”> 

where  the  a^  are  to  be  determined.  The  Schr£dinger  equation  becomes 
Ht  =  Sf  =  [flKJ  +S  )  +  IS  L  +  £  (LV  +  L~S+)]t 

Z  Z  Z  Z  m 

n 

+  2  ^ai/2(M+  2)  Y^M+l/2  X-l/2  +  “  -1/^M‘  2>  Y(4l-l/2Xl/2^ 


l  ap  Ns 


(76) 


Then 


where 


from 


‘1/2 


E  =  .1/a{P(M+  |)  +  |  (M-  ±)}  + 


■1/2 


M»-  |>  | 


*-1/2®  *  *-i/2  t«*-  j>  + 1  «<■  f  J  ♦  *1/2  »("  *  |)  | 


(77) 


B(M> 


- 1>  =/«  tj)*  -  «* 


=  B(M  *  ^) 


L*  Y  '  +  m)U  ±  m+1)  Y 


•t»+l 


Now  the  two  equations  in  a.  and  a  are  homogeneous  and  the  deter¬ 
minant  of  the  coefficients  of  a^^  must  vanish  for  a  nontrivial  solu¬ 


tion.  Hence 
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$ 


I 


& 

? 


\ 


P(M+  i)  +  |(M  -  i)  -  K 

B(il  +  |)  f 

so  that 


B(M 


5> 


x 

2 


p(M  -  |)  +  |  (M  +  |)  -  E 


=  0 


M)  =  PM  -  £  ±  ^s/p8  +2p  MX  +  X8a  +  “)*  (78) 

When  X  =  0  this  gives  the  unperturbed  energy  found  in  case  I  and  when 
P  s  0  the  unperturbed  energy  found  in  case  II. 


The  energy  can  be  expanded  for  small  X  to  compare  with  the  results 
of  perturbation  given  in  table  I .  The  expansion  to  second  order  in  X 
are 


Kd,  f) 

E(l,  £) 
Ed,  f) 
Ed,~§> 
E(l,-|) 
Ed,-f) 
E(0,  |) 
E(0,-i) 


2P  +| 


X  X? 
“  2  “  2p 

Xa 

P  +  tp 

-fi  -  ^ 

P  2p 

x  yfi 

"  2  +  2p 

-2P  +  | 


(79) 


which  are  identical  with  the  results  shown 


When  the  exact  energy  is  expanded  for 
second  order  in  p,  we  have 

E(l,±  §)  »  +  2p  +  | 

E(l,  |) 


E(1  .  I)  3  _  §  _  h  ±  1 


X  T 


§ 


E(°'  *  «)-  ±  p 


in  table  I. 

small  p  keeping  terms  through 


27  X 
_4  gf. 

27  X 


(80) 
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and  we  see  that  the  energies  agree  with  the  result  obtained  by  second 
order  perturbation  theory  as  given  in  table  II. 


The  exact  result  can  be  computed  for  a  complete  range  of  variables 
by  making  the  following  substitutions:  £  =»  x  X  ■  2(l-x)  0  <  X  <  1 
then 

E<1,  f)  =  1  +  x 

E(l,-f)  =  1  -  3x 

E<1,  |)  =  x  -  |  ±  {2xa  -  4x  + 

Edr  |)  =  “  \  ±  f3xa  -  5x  +  f)^ 

E(0,±|)  =  ±  x 

Such  a  substitution,  in  effect,  turns  off  one  of  the  interactions  as 
the  other  is  turned  on  so  that  at  the  point  x  =  0,  there  is  only  the  spin- 
orbit  interaction,  and  at  x  =  1  there  is  only  the  magnetic  interaction. 

The  proportionality  constants  are  adjusted  to  give  integer  splittings  at 
both  ends.  The  splittings  for  arbitrary  X  are  plotted  in  figure  1. 


It  is  noticed  that  the  states  which  are  degenerate  at  x  =  0  move  off 
at  different  slopes,  indicating  that  the  degeneracy  is  lifted  in  first  or¬ 
der  by  a  magnetic  field.  In  contrast,  at  x  =  1,  the  states  that  are  de¬ 
generate  move  off  with  identical  slopes,  indicating  that  second  order  per¬ 
turbation  theory  is  needed  to  lift  the  degeneracy  (curvature  is  character¬ 
ized  by  squared  or  higher  power  terms).  These  observations  are  evident  in 
the  calculations  listed  in  tables  I  and  II. 
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energy 


Figure  1, 


Relative  energy  levels 
for  the  Hamiltonian. 


as  a  function  of  the  parameter  x 


2n>  +  V< r *  +  xfLz+2Sz) +  2 n  +x )  L *S,  and  for  the  principle 
quantum  numbnr  „  ,  2.  (2„.  „  e  .  tot>l  of 

Th.  labels  1,  parenth.,..  d,„o,.  th.  orMtu  ,wUr 
•ud  tot.!  angular  momentum  projection  (M)  for  each  .tat. 
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APPENDIX  A 


Orbital  Angular  Momentum  Eigenfunctions 

* 

The  eigenfunctions  of  the  orbital  angular  momentum  operator  (squared) 
are  the  spherical  harmonics  denoted  by  Y,  (9,w) .  We  have 


L>  vt.(e,W  -  lU* i>  »,„«>, « 

(Al) 

Yt„(9’,p>  ■ ■  V9-*' 

(A2) 

where 

L8  "  '  sin8  9  W  +  ®  3®  8lR  *  9® 

(A3) 

and 

Lz  =  _i 

(A4) 

In  the  usual  manner,  the  solutions  of  La  can  be  split  into  products 
of  solutions,  one  set  being  functions  only  of  0  and  the  other  set  being 
functions  only  of  <p. 


The  functions  of  9  are  the  associated  Legendre  Polynomials  (of  the 
first  kind),  the  P^O),  and  satisfy  the  equation 

-5P5  -  5  v  -  0  <“> 

where 

p  ■  cos  9 

The  functions  of  cp  are  the  exponential  oscillatory  functions,  the 
lmcp 

♦  *  e  and  satisfy  the  equation 
da  * 

+  m8  =  0  (A6) 

where  the  m  in  this  equation  is  the  same  as  in  (AS)  and 

l  =  0,1  2,  etc.,  -  l  <  ®  <  l  (A7) 


if  the  product  solution  P.  $  are  to  have  physical  meaning. 

m 
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The  solutions  P,  (|_l)  are  given  by  a  simple  formula  in  terms  of  the 
•l® 

solutions  for  m  a  0,  the  P  (^).  These  latter  solutions  are  denoted 
simply 


V“>  ■  V“> 


(AS) 


and  are  the  so-called  Legendre  Polynomials.  They  are  polynomials  in  p 
of  degree  and  are  given  by  Rodrique's  formula 


V*  -  F7F 


d*(ua-l)^ 

.  t 

dH 


For  positive  m  the  P,  (u)  are  defined  by 

.  L/2  dm 

V*0  a  (1^  —f 4r- 

dn 


(A9) 


(A10) 


The  P^_m(li)  (m  positive  again)  are  defined  by 

■Wp)  - 


(All) 


The  spherical  harmonics  are  proportional  to  with  the  following 

normalization  (m  position  again), 


Y,  * 
tm 


(  <^-®> A,: 1/8 

\  47T  ((,+m)  If 


From  (A12)  and  (All),  one  can  deduce  that 


imcp 


(A12) 


V.- 


(A13) 


Another  useful  and  Important  property  of  the  Y.  is  their  transformation 

■o® 

under  space  inversion.  Under  such  an  inversion,  0  —  ir  -  8  and  cp  -*  cp  +  tt 
so  that 

PY^m(e,cp)  ■  Ylm(TT  -  0,  Cp  +  TT)  =  ( - )l  Y^m(®,Cp)  <A14) 

•  where  P  represents  the  inversion  operation  (called  the  parity  operator) . 
From  (A14)  it  is  clear  that  the  orbital  angular  momentum  quantum  number 
alone  determines  whether  a  state  has  pven  (positive)  or  odd  (negative) 
parity,  depending  on  whether  l  is  even  or  odd. 
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Other  useful  relations  between  the  functions  defined  above  are 


+4 


s»  -  l .  afe  A  V*  A 


■l 


*£r  V«> 


m=-^. 


where  cos  0  =  r  •  The  unit  vectors  used  as  arguments  of  the  spheri¬ 
cal  harmonics  represent  the  angle  variables  of  the  unit  vectors,  i.e., 

where  ^  a  (0,  (p) . 


We  also  have. 


IW  l 


-1—  =  Y  P.(cos  0) 
r  _tf  I  L  / +!>  I 


where  cos  8  *  ^  •  r  and  r  <  denotes  the  magnitude  of  the  lesser  of 

1^1  and  |r2|,  and  r>  denotes  the  magnitude  of  the  greater.  From 
(A15)  we  have 


47V  r<  Y.  *<r,)  Y._(r0) 


_  24.+1  t+1  ‘tm  '  1'  ^.m  2 

r,  ”  r„  .  r  > 

1  1  2  1  m 


Another  useful  formula  is 


elk‘r  =  Y  (24+1)  Ci)4  J^(kr)  V008  ®> 

A  A 

where  cos  0  m  k  •  r,  and  the  i  (kr)  are  the  spherical  Bessel  functions 
(ref  9).  Again  from(A19)  ,  ^ 

■4  | 

eik’r  =  )  4ir(i)^  j.(kr)  Y.  *(r)  Y.  (is) 

u  C  viii  -Gm 

The  formulas  given  above  are  extremely  useful  because  of  the  nor¬ 
malization  of  the  Y^m.  They  have  the  property  that, 


(A15) 


(A16) 


(A17) 


(A18) 


(A19) 
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(A20) 


where 


i&r 


^  an  *  J  d«p  J  d(cos  e> 


0  -1 

The  normalization  of  the  Legendre  polynomials  is  such  that 

.1 


I  V0  -  S^T  «u' 


A  further  property  of  the  Y,  is 

-Cni 


Vn/^  ■  I <y°>  */(p)  iL<o)> 


(A21) 


-t(m)  l'(m') 


(A22) 


where  the  C-G  coefficients  are  defined  in  Appendix  B.  From  (.A22)  and 
(A20),  one  may  obtain, 

<J  an  v^'ce.cp)  yt.(*,<p>  It ,.,<•,*)  -  ■  (?&$yi)jW Ctv»t'(0) |L(0)> 


(  L(M)> 


(A23) 


I 

f 

f 
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APPENDIX  B 


Clebschr Gordon  Coefficients 


Some  of  the  useful  properties  of  the  Clebsch-Gordan  (C-G)  coeffi¬ 
cients  are  listed  here  for  convenient  reference.  The  symbol  used  in  the 
text  is 


Some  other  symbols  for  the  same  coefficient  that  are  used  in  the  liter¬ 
ature  are 


C<W;,W>  E  6m1+m„m  C(J1J2J  iarm2  m2) 

1  « 


j,j2J 

c  1  ;  c 

"iV 


J]LJ2(Jn;mim2)  '*  <JimiJ2m2  lJmJlJ2> 


These  coefficients  are  defined  to  be  real.  Their  most  important 
property  is 

I  j(m))“  «m1+m2,m  A(Jlj2J)  (B3) 

where  implies  that  | J  7J  (  <  J  <  |da+j6)  f°r  al»y  combination 

of  a,p,Y.  This  is  called  the  triangle  condition,  nence  the  symbol 
A(J1J2J  )  . 


Other  properties  are 


<yJl<mi)J2<,n2)lJ<m>/a  (_>  <’Jl(“m,>j2("m2>lJ('m)/  (B4) 

Jl+J2-J  V 

=  <“)  |j(«")  /  (B5) 


"  (">  2  2teV8<J<-m>J2(m2)lJl(-”l)> 

CB6) 

=  (_>  1  1{S^WVJi^)J<“"),Jal“,na>^ 


(B8) 

/6__,  (B9) 

mm 
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APPENDIX  C 


The  D  Elements 

The  (r)  have  the  properties  characteristic  ol  unitary  irreducible 
representations  of  a  group,  some  of  which  are  listed  here. 


=  V  ^  ^ 


“M  ,  x 


X"S'X<'>  oJXM(r/>  -  Di'M  (rr'> 

X 

t>*T  "£<*>  DX'H'<r)  -  iSl  ««'  V  6XX' 

.  M  1  Jtrr 

where  (h  doj^  ■  I  da  I  d(cos  3)  dY  and  r  =  r  (a  3  Y) . 
o  -r 

If  r  a  r(a  3  Y)  and  r'  =  r'  (a'3  'y'),  then 

1  4<r)  4  (r'>  - 

JXM 

If  r  =  r  (o^O)  and  r  *  =  r  '  (  a'3'0)  ,  then 

OO  OO  ' 

L  47 r  dxm  (ro}  dxm  (ro')  =  6(0-n') 


One  may  extract  the  ct  and  Y  dependence  from  D  and  write 

/  X  -iM'a  /Q.  -iM'Y 

DM'M(r)  "  6  dM'M  (P)  6  (C 

where  the  d**\  are  real.  Hence,  for  a  =  Y  =  0,  the  formula  progressing 
logically  from  (C5)  and  (C6)  can  be  written 

^  ~£”  d^(3)  ^(3')  *  6  (cos  3-cos  3')  (C 
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J  J 

Other  properties  of  the  d  are  deduced  easily  from  those  of  the  D  . 

Some  of  these  are 


■4'h  <«  ■  <£>•  <»  ■  a«-  <-«  ‘ 


f* d(co«  &) 


<.«»  ■>*  -  ijfl  •„< 


(C9) 


(CIO) 


Comparing  (CIO)  with  (A21),  and  recognizing  the  reality  of  the 
d  O)  and  the  P  (p),  we  deduce  that 


P^(cos  p)  =1  d£0  O)  (Cll) 

Similarly,  comparity  (C4)  with  (A20),  it  is  at  least  consistent 

that 

“io <a  ?  °>  *  <«*> 

The  generalizations  of  (A22)  and  (A23)  are 


DX'm,(r)  Dtn(r)  "  <*'>*<*>  I L(X'+X)^>  <V (m'R(m)  |L(m'+m)^  d£,+x  m,+m  (r) 

h  (C13) 

hr  DAM(r)DXfc^E)DXm(r)a  |L<A)>  <V(m'U(m)  |L(M)> 

(C14) 
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